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Abstract 

We study the azimuthal angle decorrelation of forward jets in Deep 
Inelastic Scattering. We make predictions for this observable at HERA 
describing the high energy limit of the relevant scattering amplitudes with 
quasi-multi-Regge kinematics together with a coUinearly improved evo- 
lution kernel for multiparton emissions. 



1 Introduction 

In previous publications [U [2l [3] we have studied the azimuthal angle decor- 
relation between Mueller-Navelet jets at different hadron colliders within the 
Balitsky-Fadin-Kuraev-Lipatov (BFKL) formalism 4 beyond leading order 
accuracy [5.^ . In the present work we extend these studies to predict the decor- 
relation in azimuthal angle between the electron and a forward jet associated 
to the proton in Deep Inelastic Scattering (DIS). When the separation in ra- 
pidity space, F, between the scattered electron and the forward jet is large and 
the transverse momentum of the jet is similar to the virtuality of the photon 
resolving the hadron, then the dominant terms in the scattering amplitude are 
of the form ~ (asi^)". These terms can be resummed to all orders by means of 
the BFKL integral equation. The calculation for this process is very similar to 
that of MucUer-Navelet jets, the only difference being the substitution of one 
jet vertex by a vertex describing the coupling of the electron to the NLO BFKL 
gluon Green's function via a quark-antiquark pair. 
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This observable was previously investigated in the leading logarithmic ap- 
proximation (LO) in Ref. [6]. In the following we build on that work while still 
using the LO approximation for the jet vertex and the virtual photon impact 
factor. However, we improve the calculation by considering the BFKL kernel 
to next -to-leading (NLO) accuracy, i.e. we also keep ~ as{asY)^ terms in 
the gluon Green's function, a process-independent quantity which governs the 
dependence on Y of the cross section, previous analysis in this direction can 
be found in Ref. [7]. The NLO terms in the impact factors should be relevant 
only in the regions of moderate Y, to include them in our analysis one could 
use Ref. [5] for the leptonic part and Ref. 9J for the forward jet vertex. We will 
leave this task for future work. 



2 Forward jet cross section 

For the production of a forward jet in DIS it is necessary to extract a parton 
with a large longitudinal momentum fraction xfj from the proton. When the 
jet is characterized by a hard scale it is possible to use conventional collinear 
factorization to describe the process. Consequently, the jet production rate may 
be written as 

= J dXFj /cff(xFJ, /iF)o-(s), (1) 

with a{s) denoting the partonic cross section, and the effective parton den- 
sity [TU] being 

fes{x,fi%) = G{x,fil) + 1^ [Q^(x,/4) + Q/(x,/4)] , (2) 

/ 

where the sum runs over all quark flavors, and fj,p stands for the factorization 
scale. 

At partonic level we show a typical configuration contributing to fT(s) in 
Fig. [TJ At the leptonic vertex, we treat the quark-antiquark pair inclusively, 
while we focus on the outgoing electron, which carries momentum fci, and the 
gluon, which couples to the Green's function with momentum qi. In our notation 
the azimuthal angle of ki is ai and that of gi is 9i . 

We also work with commonly used DIS variables such as the proton mo- 
mentum P, the photon's momentum q^, its virtuality — — g^, the Bjorken 

scaling variable xbj = and the inelasticity y — p^^-^p^. Making use of 

the relation = {1 — y)Q^ and the specific structure of the leptonic vertex, we 
can write the partonic cross section in the form 

Hs) = ^ j '^'^^ J ^'1^2 J ^e"^(ki||.iepto„ic/c.<I'jct|k2> (3) 

where the rapidity is defined as y = In xpj /xbj . In bold we denote the transverse 
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Figure 1: Kinematics of the partonic cross section. 



Euclidean momenta. We can further decompose the integration and write 
ais) = ^ f; jda,jdyj<f),,j^Jd\,jd\,jdvjdv' 

n,n' — — OO 

X (y, ai |$ioptonic|qi) (qi n) {n, v\fu,\v' , n') {n' , i/'|q2) (q2 l^jct |k2)e"^. (4) 

We Irave introduced a Fourier expansion on conformal spins n, to be defined 
below. The integrals in transverse momenta are taken over the whole two di- 
mensional space while the i' integrations go from — oo to oo. The contour in the 
w-plane is to be taken to the right of all possible singularities. 

In Eq. ^ we have used the transverse momentum representation defined by 

qilqi) =qi|qi) (qi|q2) =s^'^'^ (qi - q2) , (5) 

where the kernel in operator form /C, 

/C(qi,q2) = (qi|/C|q2), (6) 
defines the BFKL integral equation at NLO, i.e., 

u;L^i+ICU. (7) 
To change representation we introduce the basis 

(qi|^,n) = -i=(q?)'''-^e»«S (8) 
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where 9i is the azimuthal angle of qi . The normaUzation in this new basis reads 



In this \n, v) representation the eigenstates of the LO kernel are 

Kq\v, n) = asXQ ^ + iv^ \v, n), (10) 
with ds = asNc/i^ and 

Xo(n, 7) = 2V'(1) - V- (7 + ^) - V- (1 - 7 + I) , (11) 

where '4^{x) = T' (x) /T{x), with T being the Euler gamma function. 

Unfortunately, to the best of our knowledge, the azimuthal angle correlation 
between the electron and a forward jet has not been extracted from the HERA 
data so far. For a future comparison with the experimental results in this work 
we implement the same kinematical cuts and constraints as those used at HERA. 
The ZEUS [TT] and HI [T^] collaborations have imposed upper and lower cuts 
in the transverse momentum of the forward jet taking into account the photon 
virtuality . They performed these cuts in order to ensure that both ends of 
the gluon ladder have a similar characteristic transverse scale. More in detail, 
they imposed 

ZEUS : i < M < 2, (12a) 

HI : ^ < ^ < ^^2b) 

These requirements are intended to suppress DGLAP evolution without affect- 
ing the BFKL dynamics. The implementation in our jet vertex of these con- 
straints is straightforward. For the ZEUS condition we have 



\j dklj d\2 (n',i.'|q2)(q2|^jct|k2 
e~™'°2 1 1 



-in a2 

Ztt 



In the case of the HI condition the 1/4 should be replaced for a 1/10. For 
simplicity, in the following, we follow the ZEUS cut. 
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In Ref. [B] it was shown that the Icptonic vertex, in our notation, reads 

cptonic I H.1 / 



,^2 



1 



d^dC 



X < I ^-e(i-0-C(i-C) + 2e(i-0C(i-C)jy' 
+ (i - 2^(1 - - 2C(i - + 12^(1 - OC(i - O) (1 - y) 



- m - 0C(1 - 0(1 - y) cos (2(01 - ai)) I 

[2aJ'"(q?,2;,Q2) + 2al''(q?,y,g2)cos (2(0i - ai)) 



(14) 



where a denotes the electromagnetic fine structure constant and J2q is the 
sum over the electric charges of the produced quark-antiquark pairs. 

To construct our cross sections we need to find the projection of this leptonic 
impact factor onto the ji', n) basis. We obtained 



d^qi (y, ai I $ioptonic | qi ) (qi 1 1^, n) 



with 



x/2 



dq? 4"\q^y,Q^)(q?)' 



(15) 
(16) 



To calculate these coefficients we need integrals of the type 

^^■ t(i-oQ' .ai-o.r 

'''£h„) ^ (^ + + "4 + + k-'-^-\ + k- ■ (17) 
with B representing the Euler beta function. Using this formula for t^/^ G {0, 1} 
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we obtained 



4°^ {^,y.Q') = 



{0" 



41/2 + 9 



tanh(7r:/) 



161/(1/2 _^ 1) cosh(7riy) 



j/2 + (i2z.2 + ii)(i_y) 



a2y2 



1 



tanh(7rz^) 



16j/(i/2 + 1) cosh(7rzy) 



X (-(4i/2 + l)(i_y)). 



(18a) 



(18b) 



The last piece needed to complete Eq. (|4]) is the gluon Green's function, 
which can be written as 



{n,v\f\v' ,n'\ 



duj 



with the eigenvalue of the BFKL kernel being 



1 



1 



(20) 



The action of the scale invariant sector of the NLO correction, given by the 
function XiI^jT); '^^^ calculated in Ref. [13]. The last term in this equation 
stems from the scale dependent part of the NLO kernel, i.e. from the running 
of the coupling. Its explicit form, in our representation, depends on the impact 
factors as given below and is discussed in more detail in Refs. [H [3]. 
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31n(2) 



2^2+1 



5 — 4cos(i/ In 4) 

cosh(27ri/) - 3 4 sin(iy In 4) ln(2) 
sinh(27rz^) 5 — 4cos(i^ In 4) 

2^6(l-y) 



V 



V + 22(1 -y)+ 4(?;2 + 6(1 - y))u'^ 
3z/2 + 1 V 



hie' y> • 



+ 1) i.2 + 1 

1 31n(2) 



(21a) 



2iy2 + i 5 - 4cos(i/ ln4) 

cosh(27ri/) - 3 4 sin(i/ In 4) ln(2) 



sinh(27rz/) 



5 - 4cos(!^ In 4) 



3zy2 



1 



3i/ 



I/(l/2 + 1) 



(21b) 
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Blending together the leptonic vertex of Eq. P^]) . the Green's function of 
Eq. (fT9)l and the forward jet vertex given in Eq. (|13p . we obtain the cross section 
of Eq. (HI , which can be expressed in differential form as 



da 



dydQ^ 



dv I dv' 



4°) {u,y,Q^) (0,H/k',0)c2(;.') 



+ {v, y, Q2) (2, j.|/>', 2)c2(i.') cos 20 



(22) 



where we have introduced the azimuthal angle (j) = a2 — oii between the electron 
and the forward jet. We have also made use of the relation {n,v\f\v\n) = 
{-n,iy\f\iy',-n). 

It is more convenient to write Eq. as 

'^'^ (y,g^r) +5(2) (y,Q2,y) cos 20] , (23) 

read 

diy (i/, y, Q2) C2(j^)e^"=x''(l"'^''), (24) 



d</> dg2 2 
where the coefficients B^") at LO read 

and at NLO: 



4lo (y^Q^Y) 



dv A(") (i/, y, ( 



X e 



s(Q')Y(xo(|«|,i^)+a,(Q=^)(xi(|nKi')-s^Xo(n,i+»'^)'irc'(!^,y.Q'))) 



(25) 



The BFKL resummation presents an instability in transverse momentum 
space when the NLO corrections are taken into account [M]. A prescription 
to increase the convergence of the perturbative expansion is to improve the 
original calculation by imposing compatibility of the scattering amplitudes with 
the collincar limit dominated by renormalization group evolution [151 116] . In 
recent publications [2J [3| we have introduced these collinear improvements to 
describe azimuthal angle dependences in the context of Mueller-Navelet jets. 
Our results were later on reproduced in Ref . . 

From a technical point of view, the collinearly-improved kernel of Ref. [5] 
differs from the one needed in the DIS case only in the term due to the running 
of the coupling in Eq. (PO)) . This contribution changes the single and double 
poles of the original kernel in the form 



ao 



02 



ao 



02 



bn 




Po /107 



(26) 

(27) 
(28) 



These equations hence replace Eqs. (24, 25) of Ref. [2]. 
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3 Phenomenology 



Besides the particular experimental cuts in the forward jet taken into account 
when calculating the jet vertex of Eq. ([13]), we also used the following experi- 
mentally motivated [18j constraints in the leptonic sector: 

20 GeV^ <Q^ < 100 GeV^ 
0.05 <y< 0.7, 

5 • 10-3 > XBj > 4 • IQ-^ (29) 
The final expression for the cross section at hadronic level reads 

da 



dY 



C^{Y)+C2{Y)cos2(i3, (30) 



with 

CniY) = 



dxFjdQ^dyU{xF},Q^)B^''\y,Q^Y)s(xFj~^ 

where we performed the convolution with the effective parton distribution of 
Eq. The index in the integral sign refers to the particular cuts of Eq. (j^ . 
The integration over the longitudinal momentum fraction xpj of the forward jet 
involves a delta function fixing the rapidity y = In xp j / xbj . It is noteworthy 
that any additional experimental upper cut on xf] would modify the coefficients 
C„, with a negligible change in their ratios. 

Since the structure of the electron vertex singles out the components with 
conformal spin and 2, the number of observables related to the azimuthal 
angle dependence is limited when compared to the Mueller-Navelet case. The 
most relevant observable is the dependence of the average < cos2(/) >= C2/C0 
with the rapidity difference between the forward jet and outgoing lepton. It is 
natural to expect that the forward jet will be more decorrelated from the leptonic 
system as the rapidity difference is larger since the phase space for further gluon 
emission opens up. This is indeed what we observe in our numerical results 
shown in Fig. [2l We find similar results to the Mueller-Navelet jets case where 
the most reliable calculation is that with a coUinearly-improved kernel. The 
main effect of the higher order corrections is to increase the azimuthal angle 
correlation for a given rapidity difference, while keeping the decrease of the 
correlation as Y grows. It is interesting to point out that, even for very small Y", 
the inclusive quark-antiquark pair (produced to couple the electron to the gluon 
evolution) generates in the case of no gluon emission some angular decorrelation 
between the forward jet and the electron. 

Finally, we estimate the theoretical uncertainties derived from not including 
the NLO impact factors by varying the scale sq, and those related to the running 
of the coupling by doing the same with the renormalization scale fi. The range 
of variation in both parameters is between 1/2 and 2 and the result is shown in 
Fig.H 
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Figure 2: < cos 20 > at the ep collider HERA at leading (solid), next to leading 
order (dashed), and for resummed kernel (dash-dotted). 




Figure 3: < cos2(^ > at the ep collider HERA for resummed kernel (dash- 
dotted). The gray band reflects the uncertainty in sq and in the renormalization 
scale II. 
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At present, the data taken at the HERA coUider provide the only possibihty 
to experimentally test our prediction. Nevertheless there are proposals to up- 
grade the Large Hadron Collider (LHC) at CERN to a Large Electron Hadron 
Collider [Hj. The projected center of mass energy y/s — lA TeV is more than 
four times bigger than at HERA and would allow for a larger rapidity separation 
between the electron and the forward jet. We use the same cuts as for HERA 
(Eq. (p9)) ) apart from an adjusted lower bound for a;Bj of 2 • 10~^. Fig. [4] is a 
plot of our results, which are very similar to those presented in Fig. [3l 



C2 
Co 




2 4 6 8 ~ 10 

Figure 4: < cos 2(/) > at a possible ep collider LHeC with a gray band reflecting 
the uncertainty in sq and in the renormalization scale ^. 



4 Conclusions 

We have studied the effect of higher order corrections to the BFKL equation on 
the angular decorrelation of forward jets in Deep Inelastic Scattering. The effect 
of these additional terms is similar to the previously studied case of Mueller- 
Navelet jets at hadron colliders. As the rapidity difference between the outgoing 
lepton and the forward jet increases, the two systems decorrelate in azimuthal 
angle due to the extra emission of soft gluons and higher order terms largely 
increase the amount of correlation when compared to the leading order calcula- 
tions. It would be very interesting to extract this dependence from the HERA 
data with forward jets and study how important BFKL effects are for this ob- 
servable. If the experience at the Tevatron is valid in this case, the BFKL 
prediction will probably lie below the data, with a gradual improvement of the 
fits as the rapidity difference increases. 
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